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1 Introduction 



The rare decays — )• /i"*"//" and B —¥ K*fi~^fi~ are sensitive probes of new particles arising 
in the extensions of the Standard Model (SM) and in particular Supersymmetry (SUSY). The 
measurements of these decays provide important constraints on the masses of new particles 
which are too heavy to be produced directly. 

At large tan/3, the SUSY contributions to the decay Bg — )• fi~^fi~ is dominated by the 
exchange of neutral Higgs bosons, and it has been emphasised in many works [l -lO that 
this decay receives large enhancement, and very restrictive constraints can be obtained on the 
supersymmetric parameters. This decay is currently searched for by three LHC experiments: 
LHCb, CMS and ATLAS, and recently LHCb collaboration reported a very strong limit on the 
branching ratio of 4.5 x 10~^ 11 which is only about 15% larger than the SM prediction. This 



improved limit further constrains the SUSY parameter space. A possible signal, although with 



a low significance is also reported by the CDF collaboration 12 

The decay B — )■ K*ix^ ix~ on the other hand provides a variety of complementary observ- 
ables as it gives access to angular distributions in addition to the differential branching fraction. 
Experimentally the exclusive B — t- K*^^ decay is easier to measure compared to the the- 
oretically cleaner inclusive mode B — )■ Xsfi^fi~ . However from a theoretical point of view in 
the exclusive mode there are large uncertainties, which come mostly from the B ^ K form 
factors. Within the QCD factorisation 13 , 14 , simplifications can be made on the form factor 



description and by looking into the rich phenomenology of the various kinematic distributions, 
observables that have smaller dependency on the form factors can be defined 15,16 . These 
observables prove to be important tools to study extensions of the SM 15 -28 and they can 



provide valuable information on different sectors of the theory. The full determination of the 
angular distributions of i? — )■ K*ii^ jjT constitutes a worthwhile challenge to the present and 
future experiments. Several angular observables have already been measured by Belle, Babar, 
CDF and LHCb. The most precise measurements come from the recent LHCb analyses with 
1 fb~^ of integrated luminosity 29 . 



In this work we study for the first time the implications of the recent measurements of 
B K*fj,'^fi~ observables on constrained SUSY scenarios and update the constraints from 
BR{B, 



30 , 31 and we 



—7- n^fi^). Our numerical analysis is performed with Superlso v3.3 
study two constrained SUSY models: CMSSM and NUHM. To give some insight on the origin 
of these constraints, a mapping of the CMSSM parameter space into the Wilson coefficients 
subspace for the most relevant operators is also provided, which displays interesting partial 
correlations and hierarchies. 

This paper is organised as follows: in section 2 we present a theoretical introduction to 
the decays Bs — ?■ n'^fJ'~ and B — t- K*fi'^fi^, provide the SM predictions and estimate the 
errors. In section 3 we summarise the experimental results and section 4 contains our numerical 
analysis of the constraints on the SUSY models that are obtained from the recent LHCb results. 
Conclusions are contained in section 5. 
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2 Observables, Inputs and Theoretical Uncertainties 



The effective Hamiltonian describing tlie b — )■ si~^i transitions lias tlie following generic struc- 
turtf|@|33|: 



n 



cff 



. „ 6 10 2 



(2.1) 



where Qi and ^2 are the scalar and pseudo-scalar operators and the primed operators are 
chirality flipped compared to the non-primed operators. Physics contributions at scales higher 
than fi are summarized in the so called Wilson coefficients cI'q,{h) (a typical choice for the 
scale is for the B decays). The Wilson coefficients include contributions from all particles 
heavier than Hb = 0{mb), in the SM these contributions are the top quark and the electroweak 
bosons and in BSM possible heavy NP particles are also summarized in the Wilson coefficients. 
The local operators Of ■* and gf ■* include the long distance contributions from scales lower than 
Hb- In the SM the primed and (pseudo-) scalar operators are either highly suppressed or absent. 
The most relevant operators for our work are 



(47r)2 
(4vr) 



(An) 



O 



10 



(Air) 



(2.2) 



gi 



'sPRh){U) 



Q2 



(4vr) 



-XsPRh){hd) 



where Pr^l = (1 ±75)/2 denote the chiral projectors and is the MS mass of the h quark. A 
full list of the operators as well as the analytical expressions for all the Wilson coefficients can 



be found in |31j . The renormalisation group equations for Ci-Cg are given in 34 , for Cq^ and 
in [l] and we calculated the running of Cg which is given in Appendix [C] 



2.1 Bs fi^fi- 

A stringent 95% C.L. limit on the untagged branching ratio BR^^tag{Bs — )• fi^fi^) < 4.5 x 10~^ 
has recently been obtained by the LHCb collaboration [llj. Taking into account the precise 
measurement 35 of the fractional decay width difference between the Bs heavy and light mass 
eigenstates, i/s = ATBji^TBj = 0.088 ± 0.014, this limit on the untagged branching ratio 
can be translated into an even stronger limit on the CP-averaged branching ratio BR{Bs — )• 
= \BR{B^s y""*'/^") + \BR{Bl /i+/i-), reading 



36 



38 



BR{Bs 12+ 12') = (1 - ys)BRunt^i.{Bs -> /i+^u") < 4.1 x 10"^ at 95% C.L. 



(2.3) 



^We neglect the doubly Cabibbo-suppressed contributions of order VubV*^. 
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Lattice QCD Group Ref. 




fs 


ETMC-11 

Fermilab-MILC-11 

HPQCD-12 


40 
41 

43 


42 


232 ± 10 MeV 
242 ± 9.5 MeV 
227 ± 10 MeV 


195 ± 12 MeV 
197 ± 9 MeV 
191 ± 9 MeV 


Our choice 


234 ± 10 MeV 


194 ± 10 MeV 



Table 1: average of lattice QCD results used in this work. 



In terms of Wilson coefficients, this average branching ratio is expressed as 31,39 : 



647r2 



* |2, 



Am?', 



m 



(2.4) 



Bs 




{Cq2 - C'nJ + 2(Cio - C[ 



m,. 



In the Standard Model, only Cio is non-vanishing and gets its largest contributions from a Z 
penguin top loop (75%), and from a charmed box diagram (24%) . With the inputs of Table § 
Cio = -4.21, from which BR{Bs l^^fi~)\sM = (3.53 ±0.38) x 10"*^. The latest experimental 
limit thus severely restraints the room for new physics, and its proximity with the 2a upper 
value calls for a discussion of the uncertainties in this SM prediction. 

The main uncertainty comes from the Bs decay constant fs^, which has recently been re- 
evaluated by independent lattice QCD groups of Table [T} Their 4.3% uncertainties agree, as 
do their results within these uncertainties, so that we have chosen an average of these three 
results in what follows. This implies a 8.7% uncertainty on the branching ratio. 

Notice this range covers the recently published result Jb^ = 225 ± 4 MeV of McNeile et 
al. 44 , whose lower value and striking precision dominates any weighted average including it, 
like the one in |45] (227 ± 4 MeV) proposed by one of the authors of 44 , or the one on [46]: 



227 ± 6 MeV. The smallness of the extrapolation error in this work raises a number of new 
questions, and we prefer to keep our naive but more conservative average. This choice has 
however little effect on the new physics applications we have in mind, as these depend mostly 
on the lower end of that range. 

Another potential source of uncertainty comes from the choice of scale at which the fine 
structure constant is used in Eq. (2.4): there is a non-negligible 4% difference between the 
running MS couplings a{mb) = 1/133 and a{mz) = 1/128. If the first choice may seem 
natural, the weak couplings involved in the top Z-penguin (or charmed box) are closer to the 
weak scale, and do not run below it, as discussed in Re f. [47| . We thus take that last value, as 
well as sin^ 9 wimz) = 0.2312 in the expression of Eq. (2.4). This may seem at odds with the 
conclusion of Ref. |47|, that choosing ce{mb) minimises the EW corrections to 5 — )■ K*i~^i~. 
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niB = 5.27950 GcV 
mK' = 0.89594 GeV 


49 
49 


ms, = 5.3663 GeV 

\VfhV,*\ = 0.0403+nnnn7 
1 ^ ts\ — 0.0007 


49 
49 




mfS(mfc) =4.19t° J« GeV 
^poie ^ -^-,2.9 ± 0.6 ± 0.9 GeV 


49 
49 


^(m,) = 1.29t° ?^ GeV 
= 0.105658 GeV 


49 
49 




as{Mz) = 0.1184 ±0.0007 

a,{nb) = 0.2161 

sin^ ewiMz) = 0.23116(13) 


49 
SI 


aemiMz) = 1/127.916 

aemimb) = 1/133 

Gp/ihcf = 1.16637(1) GeV"^ 


49 
49 




fs = 194 ± 10 MeV Table 1 
fs, = 234 ± lOMeV Table [T 


TB = 1.519 ±0.007 ps 
TB, = 1.472 ± 0.026 ps 


49 

W 




fK'A^ GeV)= 0.185 ± 0.009 GeV 
ai,j_(l GeV)= 0.10 ±0.07 
a2,±(l GeV)= 0.13 ±0.08 

yB^^K' (Q) ^ Q 4;^^ ^ Q Q4g 

As,+ (1 GeV)= 0.46 ± 0.11 GeV 


50 
51 
51 
51 
52 


fK',\\ = 0.220 ±0.005 GeV 
aij|(l GeV)= 0.10 ±0.07 
02,11 (1 GeV)= 0.09 ±0.05 
^B,^x (0) = 0.292 ±0.038 
^f""^^ (0) 0.259 ±0.036 


50 
51 
51 
51 
51 




pole 

= ml 

fJ-f — \/0.5 X fib GcV 


14 


fiQ = 2Mw 



Table 2: Input parameters used in this work. 



However, the EW corrections to Cy^g, which dominate the low qf^ region of this last process, 
are opposite to the EW corrections for Cio, which controls Bg — )■ fi^n~. Having made this 
choice for the EM-coupling, we expect EW corrections to mostly absorb the remaining scale 
dependence in aEMifJ^), leaving a small, 2% uncertainty in the branching ratio. 

The remaining theoretical uncertainties are smaller thanks to the NNLO treatment of QCD 
corrections: increasing the low scale fif, (or the matching scale fiw) by a factor of 2 induces a 
1.4% (or respectively 2%) effect. 

Finally, parametric uncertainties from the top mass (1.3%), from the Bg lifetime (1.8%) and 
from the CKM element Vts (5%), will reduce in the future. Adding all these (small) errors in 
quadrature, we thus get a Standard Model prediction assorted with an 11% uncertainty: 

BR{B, fx+fx-) = (3.53 ± 0.38) x 10"^ . (2.5) 

This value is compatible with recent SM predictions for this observable, e.g. by the CKMfitter 
group fls], for which the uncertainties are reduced by the use of other flavour information. 



2.2 B K*id+id- 

Considering the K* meson to be on-shell, the differential decay distribution of the B^ — )■ K*{^ 
K~7r~^)i^i~ decay can be written in terms of three angles 6e, 6k*, 4> and the invariant dilepton 
mass squared (g^) [15 



22 



9 
32^ 



J(g^, 9i, 9k*, 0) dq'^ dcos 9i d cos 9k* 



(2.6) 
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In the above equation 9i is the angle between l~ and in the rest frame of the dilepton, 
Ok* is the angle between K~ and in the K* rest frame and is the angle between the 
normals of the plane and the K~tt^ in the B^ rest frame. The angular dependence of 
J(g^, Oi, Ok*, 0) can be written as 



J(g^, On, Ok* , 0) = JI sin^ Ok* + J{ cos^ Ok* + ( ^2 sin^ Ok* + ^2 cos^ ^i^*) cos 20i> 

+ J3 sin^ Ok* sin^ 6*^ cos 20 + J4 sin 26^^^* sin 20 1 cos + J5 sin 26'x* sin Of, cos 

+ ( Jg sin^ + Jg cos^ Ok*) cos + J7 sin 26'x* sin sin 

+ Js sin 20k* sin 2^£ sin + Jg sin^ Ok* sin^ 6'^ sin 20 . (2-7) 

The angular coefficients Jj (see Appendix jA[ ) , are functions of and can be described in 
terms of eight transversity amplitudes , Ay ' , A^ , At, As- The transversity amplitudes up 
to corrections of 0{as) can be expressed in terms of seven independent form factors, ^0,1,2, 7^1,2, 3 
and V . Since these form factors are hadronic quantities they require non-perturbative calcula- 
tions and hence are a main source of uncertainty in the exclusive mode. But even if the form 
factors were known precisely the B^ — )■ K*i~^£~ decay would still need corrections emerging 
from non-factorisable effects which are related to the current- current operators Oi and O2, the 
QCD penguin operators O3-OQ and the chromomagnetic operator Og- The non-factorisable 
corrections contribute to the decay amplitude by producing a virtual photon which decays into 
a lepton pair. When the K* has a large energy (g^ is small) and the decaying hadron (i?°) 
is heavy, the non-factorisable corrections can be computed in the QCD factorisation frame- 
work 13 14 . In the QCDf framework in the large recoil limit the seven independent form 

While a, corrections to the form 



53,54 



factors can be expanded in ratios of l/rnt and 1 / Ek* 
factors in QCDf have been calculated [54], the l/m?, corrections are unknown. Neglecting l/rub 
corrections, the transversity amplitudes at NLO in as in the large recoil limit arfj^ [l5[|22[|23| : 



[(Cg + C',) T (ClO + C[,)] 



2mb 



Arrib Ek 



Mb + rriK* <f 



(2.8a) 



[(Cg - C',) T (Clo - C[,)] 



Mb — rriK* 



+ 



Mb 



-T; 



(2.8b) 



^For At and Ag, the correction which emerges from describing Aq in terms of is very small and has 
been neglected. 
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L.R 



N 



:|[(C9-C^)T(Cio-C;o 



)] 



{Ml - ml, - q^){MB + mK*)Ai{q^) - X 



A2{q' 



Mb + rriK* 



A, 



+ 2mb 



Mb 
2(Cio 



{M'b + 3ml, - r)Tj_" 



Ml - ml. 



C[o) + 



q 



memb 



(Cq2 ~ C'q 



mx* 



^s = -'^v^X{Cc,,-C',J^^,{q^), 
mj) mx* 

where Ek* is the energy of the final vector meson in the B rest frame 



(2.8c) 
(2.8d) 
(2.8e) 



Ek* 



Ml + ml, - g2 



2M. 



B 



(2.9) 



and 



N 



r^2 2 

^ F em 

3 ■ 2i0 7r5M| 



\Vt,V:fq^^l3, 



1/2 



Amj 



(2.10) 



and A is the triangle function A(M^, m|-, , g^) 

\ = Ml + ml,+q^- 2{Mlml, + ml,q^ + Mlq^) 



(2.11) 



B 



Further explanation on what we use for T^y, Ai^2 and V is given in Appendix 

For the differential distribution of the CP conjugate decay — )■ K*{-^ K^7T^)t^i^ , if we 
keep the definition of to remain the same and replace with for the definition of 6k* 
and consider to be the angle between the normals of the plane and the K'^n' in the 
B^ rest frame, we can use Eq. (2.6) where J is replaced with J. The function J{q'^,6£,6K*,(p) 



is obtained from Eq. (2.7) by the replacements 55 



t(c,s) , j{c,s) 
■^1,2 '^1,2 ) 



t(c,s) . _ t(c,s) 
Jc r J p, , 



■^3,4,7 J' 



3,4,7, 



J' 



5,8,9 



^ -J' 



5,8,9, 



(2.12) 
(2.13) 



where J is equal to J with all the weak phases conjugated. The change of sign for J5,6,8,9 can be 
understood by considering that the CP conjugate decay leads to the transformations 9i ^ n — 9i 
and -0 in Eq. ^7f\. 
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2.2.1 Observables 



Integrating Eq. (2.6) over all angles, the dilepton mass distribution can be written in terms of 



the angular coefficients 22 , 23 



dT 



4 



Ji 



(2.14) 



For the normalised forward-backward asymmetry ApB use 22[23 

fO 







d cos 01 



dT 



dq^ d cos 9i / dq^ 




(2.15) 



The longitudinal polarisation fraction can be constructed as the ratios of the transversity 
amplitudes and therefore contains less theoretical uncertainty from the form factors. For the 



longitudinal polarisation, Fl we have 15 , 22 



Fds) 



Of) 



dV/dq^ 



(2.16) 



where Ji = 2J^' + for i = 1,2, 6. 

An interesting observable is the zero-crossing of the forward-backward asymmetry (gg). 
Neglecting the lepton masses and not considering the chirality flipped and scalar operators for 
the zero of Aps we obtain 



2 ~^br. 



(2.17) 



To calculate the zero of Aps in our numerical analysis we have directly used the relation for 



Afb (Eq. (2.15)). 



Another observable which is rather independent of hadronic input parameters is the isospin 
asymmetry which arises from non-factorizable effects. The non-factorizable effects depend on 
the charge of the spectator quark, and hence depending on whether the decaying B meson is 
charged or neutral, there will be a difference in the contribution of these effects to the decay 
width which can cause an isospin asymmetry. The (CP averaged) isospin asymmetry is defined 



as 20 



dAi 
dq^ 



(2.18) 



dT[B^ K*^£+r]/dq^ - dT[B^ K*^£+r]/dq^ 
dr[BO K*H+l-]/dq^ + dT[B^ -> K*^£+£']/dq^ ' 

the smallness of the isospin asymmetry makes it 



14 



In the SM, dAj/dq^ is less than 1% 
sensitive to isospin-violating NP effects. 

Other observables of interest are the transverse amplitudes which have a small dependence 
on the form factors and a large sensitivity to right-handed currents via C^. They are defined 



as 15,16 22 56 
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|(2)' ^ J- 



A'4\s) = , (2.19) 



2J, 



^i^^._(mi±^\" .2 21) 



AUs) - ^ , (2.22) 

which is sensitive to the complex phases, but very small ((9(10"'^)) in the SM. 

All the observables can also be expressed in terms of the CP averaged angular coefficients 



as well as 



Si 22 . In particular, 5*3, which has recently been measured by LHCb, can be related to Fl 



and Arj.' by 



S, = 1(1-Fl)AP . (2.23) 



For the high region we follow |26|, where the heavy quark effective theory (HQET) of [57j 
is used. In that framework improved Isgur-Wise form factor relations are used to extrapolate 
the lacking high form factors from the calculations of the low form factors and an operator 
product expansion in powers of 1/Q^, {Q = {mi,, is used for the estimation of the long- 

distance effects from quark loops. For the low recoil region (high g^), the following observables 



are proposed by 26 



4^) = , ^-^^ , (2.24) 
v/-J|(2J|- Js) 

4) = , ^"^^ , (2.25) 
v/-2 J|(2J| + J3) 

^ ^^'^^ . (2.26) 
2v/(2J|)2 - (J3)2 

These transversity observables have small hadronic uncertainties at high g^. 



2.2.2 Standard Model values and theoretical uncertainties 

For our numerical analysis we have used the NLO relations for the transversity amplitudes in 
the (1 < < 6 GeV^) region where we are sufficiently below the charm resonance threshold 
(g^ = m^) and far enough from the kinematic minimum where the decay amplitude is dominated 
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CM 










Ceifib) 


Cfifib) 


Cfifib) 


CgifJ'b) 


CioiiJb) 


-0.2610 


1.0076 


-0.0052 


-0.0795 


0.0004 


0.0009 


-0.2974 


-0.1614 


4.2297 


-4.2068 



Table 3: SM Wilson coefficients at fXh = ^ and /xq = 2Miy to NNLO accuracy in a, 



Observable 


SM prediction 


(FF) 


(SL) 


(QM) 


(CKM) 


(Scale) 


10^ X BR{B ir>+/i")[i,6] 


2.32 


±1.34 


±0.04 


+U.04 
-0.03 


+0.08 
-0.13 


+0.09 
-0.05 


{AFB{B^K*ij+fi-))[,^e] 


-0.06 


±0.04 


±0.02 


±0.01 






(Fi(S^ir>+/i"))[i,6] 


0.71 


±0.13 


±0.01 


±0.01 






q^{B A^>+/i-)/GeV" 


4.26 


±0.30 


±0.15 


+0.14 
-0.04 




+0.02 
-0.04 



Table 4: SM predictions and errors. 



by the photon pole. The input parameters can be found in Table |2] and the SM Wilson 
coefficients are given in Table p| The available experimental values are given for bins which 



can be shown as (observable) ,q2^^j, in other words the dq^ integration is over the [^mm' Qmaxl 
bin. For normalised quantities like ApB and Fl the numerator and denominator are integrated 



separately 15 , 26 



r,2 ^2 



O /J \ 

(^^i^)[.L„,e.J - -8(rfr/rfg2) ' (2-27) 

/ TP \ ^ '^'^ Vlmin-'^rnax] /r) r\Q\ 

We estimate the theoretical uncertainties for the SM values in the low region using two 



methods: one follows the approach of 58 , and the second is based on a Monte Carlo method. 
Both methods give similar results. We consider five different sources of uncertainty. The errors 
from the form factors (FF) have been calculated by adding in quadrature the uncertainties 
due to V, Ai and A2 (11%, 13% and 14%, respectively (5l]). For the unknown l/m^ sub- 
leading corrections (SL), we have assumed 7% corrections to the 71, y terms in the transversity 
amplitude^ these corrections have been added in quadrature. Another group of errors is 

pole 



QM) which we have estimated by separately varying m 



from the quark mass uncertainties 

and m^^ according to Table 2j and added in quadrature. Another source of error comes 
from the uncertainty in the CKM matrix element combination (CKM) which has been 

computed by considering the uncertainty given in Table |2} The last source of error that we 
consider is the scale dependence (Scale) which we estimate by varying fib between /ib/2 and 2fih 
(with /ib = m^"'*^). These five groups of errors for BR, Aps, Fl and gg in the SM have been 
gathered in Table |4j For the Standard Model predictions the primed coefficients (CIqj^^qJ as 
well as the scalar coefficients (Cqj j) have been put to zero. 



■^If for the (SL) error a 10% correction to the T±.\\ is considered the overall uncertainties will not have a 
significant change since they are dominated by the (FF) uncertainties. 
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Observable 


SM prediction 


Experiment 


lO^GeV X {dBR/dq^ {B ii'>V))[i,6] 


0.47 ±0.27 


0.42 ±0.04 ±0.04 


(A^B(5^ii'V/^"))[l,6] 


-0.06 ±0.05 


n 1 0+0.U6+U.U1 

"J- -"-0-0.06-0.01 


(FL(5^i^>V-))[i,6] 


0.71 ±0.13 


n ««+0.U6+0.U4 
'-'■^"-0.06-0.03 




^•^"-0.34 





Table 5: Experimental values and SM predictions (the theoretical errors are from adding in 
quadrature the different errors in Table |4]) . 



For the high q region we have used the relative errors of Table 2 in 26 . 



3 Experimental results 

At present, the best upper limit for BR(i?s — t- measured in a single experiment comes 

(3.29) 



fromLHCb 11 



BR{B, fi+fi-) < 4.5 X 10" 



at 95% CL. This upper limit is followed by the result from CMS 59 

BR{Bs i^^fi') < 7.7 X 10"^ . 

The CDF collaboration obtains a 95% C.L. upper limit fT2^: 

BR{Bs fi^fi') < 3.4 X 10"^ , 

together with a one sigma interval 

BR{Bs = (1.3^[J:?) X 10~* , 



(3.30) 



(3.31) 



(3.32) 



coming from an observed excess over the expected background. 

For our numerical evaluations, we consider the LHCb limit, and accounting for 11% theo- 



retical uncertainty (as explained in section 2.1), we impose the following limit at 95% C.L. 

BR{Bs fi^fi') < 5.0 X 10"^ . (3.33) 
For B — )■ K*^^ related observables we also use the latest LHCb results which correspond 



to an integrated luminosity of 1 fb~^ 29 . The results are summarised in Table [s] where the 
experimental uncertainties are statistical and systematic. For comparison, the SM predictions 
with the corresponding theoretical errors (from the five sources of errors mentioned in Table |4| 
added in quadrature) are also provided. 

In addition to the observables in Table [5} three other observables have also been measured 
using 1 fb~^ of LHCb data namely 5*3 and Aim |29] and also very recently the isospin asymmetry 
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Aj 60 . The reported results are: 

(2 53(i?^ir>V))[i,6] = O.mt'oilTol, (3.34) 
(A,„(5^K>V))[i,6] = 0.07^°:?™, (3.35) 
(A/(5 -> K*/i+/i"))[i,6] = -0.15 ±0.16. (3.36) 

However, with the current experimental accuracy, these observables are not sensitive enough 
to probe SUSY parameters. 

4 Constraints on SUSY 

To illustrate the impact of the recent limit on BR(i?5 — )■ /U+yU~) and the measurements of 
the angular distributions and differential branching fraction of the B — )■ K* ^'^ decay by 
LHCb, we consider constrained MSSM scenarios. The observables are calculated as described 



in section 2 using Superlso v3.3 30 31 



We focus on two specific scenarios, both assuming SUSY breaking mediated by gravity and 
invoking unification boundary conditions at a very high scale mouT where the universal mass 
parameters are specified. The first model is the CMSSM, characterised by the set of parameters 
{mo, mi/2, v4o, tan/3, sgn(/i)}. Here ttiq is the universal mass of the scalars, mi/2 the universal 
gaugino mass, Aq the universal trilinear coupling, and tan (5 the ratio of the vacuum expectation 
values of the Higgs doublets. The second model we consider involves non-universal Higgs masses 
(NUHM). This model generalises the CMSSM, allowing for the GUT scale mass parameters of 
the Higgs doublets to have values different from itlq. The two additional parameters can be 
traded for two other parameters at a lower scale, conveniently the parameter and the mass 
Ma of the CP-odd Higgs boson implying that the charged Higgs boson mass can be treated 
essentially as a free parameter, contrary to the CMSSM. 

First we study the constraints from BR(i?s — )■ /^^A*^) in the CMSSM. For this purpose, 
we scanned over mo G [100, 3000] GeV, mi/2 G [100, 3000] GeV, Aq e [-5000, 2000] GeV and 
tan/3 G [1, 60], with yU > 0. We generated about 500,000 points, and for each point we calculate 
the spectrum of SUSY particle masses and couplings using SOFTSUSY 3.2.4 IM^ and compute 
BR(i?5 — > yU^/i~) using Superlso v3.3. The constraints are shown in Fig. IT] in the planes 
(M(-j, tan/3) and (M//±, tan/3). These results could be compared to the previous results using 



the LHCb 2010 collected data [10,62 . It is clear that the large values of tan/3 are strongly 
constrained. 

In the NUHM scenario, the parameters are scanned over in a similar way as in the CMSSM, 
and in addition /i G [-2000,2000] GeV and Ma G [50,1100] GeV. The results are shown in 
Fig. |2} Since there are two additional degrees of freedom in NUHM as compared to CMSSM, 
it is easier for a model point to escape the constraint, as can be seen by comparing Figs, [l] 
and [2] where the allowed points are displayed on top. On the other hand since Mj^ and Mh± 
are discorrelated in NUHM, it is possible for any Mi^ to have a very small Mh± and therefore 
being excluded, as can be seen from the plot with the allowed points in the background. 
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Figure 1: Constraint from BR(i?s — )■ in the CMSSM plane (Mf^, tan /3) in the upper 

panel and {Mh± , tan (3) in the lower panel, with the allowed points displayed in the foreground 
in the left and in the background in the right. 



The effect of an SM-like measurement of BK{Bs fi'^fi ) in non-constrained MSSM (the 



pMSSM) is demonstrated in 63 



Next we consider the constraints from B — >■ K*fi^fi~ observables. In order to study the 
maximal effects we consider tan /3=50 but show also the results for tan /3=30, and investigate 
the SUSY spread in function of the lightest stop mass. We start with the averaged differential 
branching ratio as defined in Table |4j The results in CMSSM are displayed in Fig. [3] for the low 
region and in Fig. |4] for the high region, where the solid red lines correspond to the LHCb 
central value, while the dashed and dotted lines represent the 1 and 2cr bounds respectively, 
including both theoretical and experimental errors (added in quadrature). At low g^, this 
branching ratio excludes M^^ below ~ 250 GeV for tan /3=50 and ~ 150 GeV for tan /3=30. In 
the high region the branching ratio is doing slightly better, as the M^^ below ~ 300 GeV and 
~ 200 GeV are excluded for tan/3=50 and tan/3=30 respectively. As this light stop region is 
already excluded by the direct SUSY searches for the same scenario, BR{B — )■ K*fi'^fi~) does 
not provide additional information. The main reason of the limited constraining power of the 
branching ratio is the large theoretical uncertainties (mainly due to form factors) from which 
this observable is suffering. The results are shown for two values of Aq (=0 and -1000 GeV) 
for comparison. As can be seen from the figures, the constraints with smaller Aq are slightly 
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Figure 2: Constraint from BR(i?s — )■ fi~^fJ^~) in the NUHM plane (Mf^, tan/3) in the upper 
panel and {Mh± , tan (3) in the lower panel, with the allowed points displayed in the foreground 
in the left and in the background in the right. 



stronger. 

Contrary to the branching ratio, angular distributions, in which the theoretical uncertainties 
are reduced, can in principle provide more robust constraints on the SUSY parameter space. In 
particular, we consider in the following the forward-backward asymmetry ApB, the zero-crossing 
ql of ApB, Fl, as well as 5*3, Ajrn for which the LHCb results with 1 fb~^ of data are available. 
The two latter observables do not provide any constraint with the current results and accuracy. 
The SUSY spread in function of the stop mass of ApB, Qq and Fi is given in Figs. [5]-[7j As can 
be seen, ApB provides the most stringent constraints among these observables, and excludes 
~ 800 GeV and 600 GeV at tan/3=50 and tan/3=30 respectively. Qq on the other hand 
excludes M^^ < 550 GeV (for tan (3=50) and 400 GeV (for tan (3=30) while Fl excludes M^^ < 
200 GeV (for tan/3=50) and 150 GeV (for tan (3=30). The impressive constraining power of 
ApB mainly due to the fact that the measured central value is smaller than the SM prediction 
and in addition the reported experimental errors are more than twice smaller than the previous 
results. 

Same observables at high have less impact on the SUSY parameters and therefore their 
results are not reproduced here. 
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Figure 3: SUSY spread of the averaged BR{B — )• K* jj,^ n ) at low in function of the hghtest 
stop mass, for tan/3=50 (upper panel) and tan/3=30 (lower panel), in the left for = and 
in the right for Aq = -1000 GeV. 
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Figure 4: SUSY spread of the averaged BR{B — )■ K*^^fi ) at high g^, similar to Fig. pi 
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Figure 5: SUSY spread of the Afb{B — )■ K*fi^^ ) at low in function of the hghtest stop 
mass, for tan /3=50 (upper panel) and tan /3=30 (lower panel), in the left for Aq = and in the 
right for Aq = -1000 GeV. 
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Figure 6: SUSY spread of the Apb{B K*fi^^ ) zero-crossing, similar to Fig. [sj 
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Figure 7: SUSY spread of the Fl{B — )■ K*fi~^fi ) at low in function of the hghtest stop mass, 
for tan /3=50 (upper panel) and tan /3=30 (lower panel), in the left for Aq = and in the right 
for Ao = -1000 GeV. 
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Figure 8: SUSY spread of the Afb{B — )■ K*fi^fi ) at low in function of the pseudo scalar 
Higgs mass (left) and the gluino mass (right), for tan/3=50 and Aq = 0. 



To illustrate the sensitivity of the observables to other SUSY parameters, in Fig. |8] we 
present the variation of ApB (which has the largest impact as we have shown) with respect to 
the pseudo scalar Higgs mass and the gluino mass. 
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Figure 9: SUSY spread of the isospin asymmetry {B — t- K*fi'^fi~) at low in function of the 
hghtest stop mass, for tan/3=50 and Aq = 0. 
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for tan /3=50 and Aq = 0. 



Another observable of interest for which LHCb has announced the measurement very re- 
cently 60 is the isospin asymmetry Aj, defined in Eq. (2.18). Fig. |9| illustrates the SUSY 



spread of Aj in function of M^-^. Since the experimental measurement (Eq. (3.36)) has an error 
larger than the SUSY spread, the current result does not provide any information on the SUSY 
parameters. 



(2) .(3) .(4) 



In Fig. 10, we show the SUSY spread in function of the lightest stop mass, for Aj, , Aj, , A. 



(defined in Eqs. (2.19)-(2.21)) which are not measured yet by LHCb. Large deviations from the 
Standard Model can be seen for small values of M^^, and depending on their eventual upcoming 
experimental measurements they could be of interest in constraining SUSY parameters. 
A comparison between these observables in the plane 
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As 



[mi/2,'mo) is given in Fig. 
expected, Aps is the most constraining observable also in this plane. All the observables show 
more sensitivity at larger tan/3 and smaller Aq. 

To trace the origin of these constraints, it is interesting to present them in terms of Wil- 
son coefficients. The variation of the observables in SUSY is driven on the one hand by the 
additional contributions to the Wilson Coefficients, in particular Cy, Cg, Cg and Cio, and on 
the other hand by new contributions through the scalar and pseudo scalar coefficients Cq^ and 
Cq2- The full variation of the Wilson coefficients in a full CMSSM scan is presented in Fig. [T2 
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Figure 11: Constraining power of the B — )■ K*fi~^fi observables in the plane (mi/2,wio), for 
tan/3=50 (upper panel) and tan/3=30 (lower panel), in the left for Aq = and in the right for 
Ao = -1000 GeV. 



ignoring existing constraints on SUSY parameters or Wilson coefficients. As can be seen, C7 
and Cg can have both signs in SUSY and their correlation is visible in the figure. As for Cg, 
it varies only by very little while Ciq can have a larger spread. This feature can be under- 
stood once we notice that box diagrams are suppressed with respect to scalar- or Z-penguin 
diagrams, giving 6Cg/6CiQ ~ {gv/gA)fj. ~ 1 — 4sin^6'vK- On the other hand, scalar coefficients 
can receive very large contributions in SUSY, as already known. An interesting information 
can be obtained from the zero-crossing of ApB since it depends on the relative sign of C7 and 



Cg as can be seen in Eq. (2.17), and excludes the positive sign of C7 as can be seen in Fig. 13 
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64 , 65 while the latest 



Finally, in Fig. 14 we show the correlations of BR(i? — )■ Xs'y) with BR(i? — )■ K* jji^jj,' 
and Afb{B — )■ K*fi'^fj,~). In the SM, BR{B — )■ Xg'y) is dominated by contributions from C7. 
The SM prediction for this branching ratio is (3.08 ± 0.22) x 10~^ 
combined experimental value from HFAG is (3.55 ±0.24 ±0.09) x 10" 
are strong correlations between BR(i? — > ^^7) and B — )■ K*fi~^fi~ branching ratio and forward 
backward asymmetry, especially for small and intermediate values of tan/3, where the scalar 
and pseudo scalar contributions have a limited effect. 



66 . As expected, there 
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Figure 13: Variation of Cj and Cg in the CMSSM with all the parameters varied in the ranges 
given in the text. Only the green points are allowed by the zero-crossing oi ApB- The red lines 
correspond to the SM predictions. 
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Figure 14: The correlation between BR(i? — )■ ^^7) and averaged BR(i? — )■ K*fi'^^ ) in the 
left, and between BR(i? — )■ Xs7) and Afb{B — )■ K*^'^fi~) in the right. 

5 Conclusions 

The rare decays Bs — ?■ ^Ji^ and 5 — )■ K*fi~^fi~ are of great importance in exploring the foot- 
prints of physics beyond the SM. We have presented in this paper an update of the constraints 
by the recent LHCb limit on BR(i?s — )■ in two constrained MSSM scenarios, CMSSM 

and NUHM. 

We also demonstrated that B — t- K*fi^fi'' provides a large number of complementary ob- 
servables from the full angular distribution. In particular while Bs — ^ /i^yU" lose very rapidly its 
sensitivity by lowering down tan/3, the constraints by B — )■ K*n'^fi^ observables weaken only 
mildly and thus act as additional tools to explore SUSY parameters. The forward backward 
asymmetry is in particular a very powerful observable which provides competitive information 
to the direct SUSY searches. To display the origin of these constraints in the CMSSM, we 
have shown the possible values and correlations between the Wilson coefficients of the various 
operators relevant for these processes. Among those, the Wilson coefficient Cg is a parameter 
for which these angular observables provide a unique handle. However, we have shown that the 
sensitivity of Cg to CMSSM effects is suppressed by the muon vector coupling, because they 
proceed mainly through a Z penguin diagram. The constraining power of other observables, 
like the isospin asymmetry whose measurement is under way, have been shown to display their 
interest in the framework of Minimal Flavour Violation models for which the CMSSM provides 
an interesting prototype. 

With reduced theoretical and experimental errors for these observables, in particular for 
the forward backward asymmetry and its zero-crossing which have been recently measured for 
the first time by the LHCb experiment, one would have access to very powerful observables in 
constraining supersymmetric models. The information obtained in this way can also serve as 
consistency checks with the results of the direct searches. 
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A Angular coefficients Ji 



The angular coefficients are given below 22 
(2 + 



r 



J3 



4m« 



^ir + i4r + ^ ^)] + (AiAf + ) 



4m« 



[|A,p + 2Re(A^O]+/3||A5p, 
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\L\2 \ aL\2 



\A\\\' + {L^R 
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J5 = 



me 



Re{A^A*s + AfA*s[ 



Jl = 2/3, [Re(4Ai*) - (L i?)] 



Jg^ = 4/3,^ Re [AqM^ + (L -> i? 



J7 = V2 A 

Js = i=/3| [lm(A^Ai*) + (L ^ i?)] , 
Jg = [lm{AfAi) + {L^R 



Im(A^4*) -{L^R) + ^ MAiA*s + AfA*s) 



with Pe defined in Eq. rt2.10) 



(1.37a) 

(1.37b) 

(1.37c) 
(1.37d) 
(1.37e) 

(1.37f) 

(1.37g) 

(1.37h) 
(1.37i) 

(l-37j) 

(1.37k) 
(1.371) 



B Form factors 

The 71,11 (which include non-factorisable corrections) are described in terms of the soft form 
factors and can be found in [13 14 where, in order to get T^ip Cf^ appearing in 7j_,|| should be 
replaced with {Cf^ ± Cy) [23]. To obtain the soft form factors we have used the factorisation 
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V 

A2 



ri 



r2 ra 



0.923 -0.511 5.322 

0.290 
-0.084 0.342 



49.40 
40.38 
52.00 



Table 6: Fit parameters describing the dependence of the form factors V and Ai 2 in the 



LCSR approach 51 . 



scheme used in 14 : 



Mb 2^ 



Mb + rriK* 



2E 



K* 



M 



B 



(2.38) 
(2.39) 



The full form factors V and Ai 2 have been taken from light-cone sum rule (LCSR) calculations 
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A2{q 
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"fit 
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l-q^/m)^^ ' 
2^ ^1 
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^2 



1 - q^/m)^, (1 - gV"^/iJ^ 



(2.40) 
(2.41) 
(2.42) 



where the fit parameters ri 2,m'j^ and rnj-^ are given in Table [6j 

C Renormalisation group equations for Cg 

To obtain the evolved Wilson coefficients from the matching scale (/xq) to Hb, the renormalisation 
group equation (RGE) for C has to be solved: 



(3.43) 



where 7 is the Anomalous Dimension Matrix (ADM). A general solution for the RGE is 

C(/ife) = ?7(/i„,/io)C(/io) , (3.44) 

where U is the evolution matrix 

f/(/i,,^o) = T.exp / dg'^-]-}^ . (3.45) 
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Tg is the g(= A/47ra^) -ordering operator such that the couphng constants increase from right 
to left. 7 and /3 can be expanded in g 
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The Wilson coefficients and the evolution matrix expanded in are 



f/(/ib,/io) = I 



An 



^^^^ Q (^). 



47r 



Following the conventions used in |34j , the evolution matrix is 

n 9 



in) 
kl 



j=0 1=1 



where rj = a^(/io)/as(Aife)- 

Using the above relations for the different orders of Cg{fii,) we have 



=1 1=1 



.«=i i=i 
To obtain the evolution matrix we have followed 



1=1 



(3.46) 
(3.47) 



(3.48) 
(3.49) 



(3.50) 
(3.51) 



(3.52) 
(3.53) 
(3.54) 
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and 68 . Taking 7^°), 7(1) and 7(2) 



from 34 and [69], we have produced the necessary "magic numbers" {mku) for the evaluation 
of f/g;. Since O9 does not mix with Oy g the magic numbers mgyj and mggj are all zero. The a^'s 
and the rest of the mku can be found in Tables [7p!3 
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where to 



In the above formulas, Ci"'\ijq) for Z 7^ 7, 8 can all be found in section 2 of 
get our c[^\^o) we use C[^\^q) = Cf^\^Q) — Cf"'\fio), and C'^l{^o) are not needed since t/gr 
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Table 7: The numbers aj used in Eq. (3.50) 
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Table 8: "Magic numbers" Tn^i^\ 



and Ugs are zero. To get the Wilson coefficients based on our choice of operators (Eq. (2.2)) 
we have 
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Table 9: "Magic numbers" Tn^i^\ 
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^93i 














-0.7405 


-0.1088 


0.3825 


-7.9139 





(11) 














-0.4942 


0.0856 


0.3111 


0.6465 





(11) 














-14.3464 


-2.2495 


1.8402 


-53.6643 





(11) 

^9ei 














-6.6029 


1.0948 


0.0584 


-3.2339 





(11) 
^99/ 






























Table 10: "Magic numbers" m. 



i 


1 


2 


3 


4 


5 


6 


7 


8 


9 


(20) 








0.6878 


-0.9481 


-0.1928 


-0.8077 


-0.2554 


0.0562 


-0.6436 


(20) 








1.3210 


3.1616 


-0.4814 


1.9362 


-5.0873 


0.0468 


-13.5825 


(20) 
^93i 














-2.5758 


-5.8751 


0.0922 


0.6433 


7.7756 


(20) 














-2.6194 


1.1302 


-27.7073 


-0.8550 


16.0333 


(20) 
^95i 














-6.4519 


-555.931 


35.1531 


80.2925 


102.043 


(20) 














-53.3822 


34.3969 


-124.609 


-32.7515 


-98.8845 


(20) 
^99/ 






























Table 11: "Magic numbers" m. 



i 


1 


2 


3 


4 


5 


6 


7 


8 


9 


(21) 
"^91/ 








-1.7394 


0.0530 


0.1741 


-0.1036 


-0.1478 


1.2522 





(21) 








-2.5918 


-0.2971 


-0.5949 


0.1241 


0.3170 


2.8655 





(21) 
^93/ 














1.6188 


-0.5311 


-10.9454 


4.36311 





(21) 
"^94i 














0.1967 


-0.3775 


0.7915 


15.2328 





(21) 
^95/ 














-74.1049 


-32.6399 


-35.4688 


607.188 





(21) 
^96i 














29.5971 


-1.7519 


50.6366 


216.094 





(21) 






























Table 12: "Magic numbers" m. 



26 



i 


1 


2 


3 


4 


5 


6 


7 


8 


9 


(22) 
^91i 








4.1531 


-0.4627 


-0.3404 


-1.0326 


0.0809 


0.2167 





(22) 
"^92i 








12.4592 


0.6940 


-1.7340 


1.2360 


-0.1921 


0.3998 





(22) 














-7.6198 


-4.1683 


1.5484 


15.674 





(22) 
^94i 














-5.0848 


3.2810 


1.2592 


-1.2804 





(22) 














-147.615 


-86.199 


7.4486 


106.285 





(22) 
^96i 














-67.9394 


41.9523 


0.2364 


6.405 





(22) 
"^99/ 






























Table 13: "Magic numbers" m. 
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